Fermion mediated long-range interactions of bosons in the ID Bose-Fermi-Hubbard model 
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The ground-state phase diagram of mixtures of spin polarized fermions and bosons in a ID periodic lattice is 
discussed in the limit of large fermion hopping and half filling of the fermions. Numerical simulations performed 
with the density matrix renormalization group (DMRG) show besides bosonic Mott insulating (MI), superfluid 
(SF), and charge density-wave phases (CDW) a novel phase with spatial separation of MI and CDW regions. We 
derive an effective bosonic theory which allows for a complete understanding and quantitative prediction of the 
bosonic phase diagram. In particular the origin of CDW phase and the MI-CDW phase separation is revealed as 
the interplay between fermion-induced mean-field potential and long range interaction with alternating sign. 

PACS numbers: 



Ultra-cold atomic gases in light induced periodic potentials 
have become an important experimental testing ground for 
concepts of many-body physics since they allow the realiza- 
tion of precisely controllable model Hamiltonians with widely 
tunable parameters. A system which has attracted particular 
interest in the recent past is a mixture of bosons and spin- 
polarized fermions in a deep lattice potential, described by the 
Bose-Fermi-Hubbard model (BFHM) [1, 2]. Mixing lattice 
bosons with fermions, Giinther et al. [3] and Ospelkaus et al. 
[4] observed an unexpected reduction of bosonic superfluid- 
ity which triggered a number of theoretical and experimental 
studies on the influence of fermions on boson superfluidity. In 
the limit of small fermion mobility the phase diagram can be 
well understood by mapping to the purely bosonic Hubbard 
system with binary disorder [5, 6]. For increasing fermionic 
hopping amplitudes a number of new phenomena emerge [7], 
including polaronic phases [8] and density waves [9, 10]. Fur- 
thermore under certain conditions bosons can enter the a su- 
persolid phase (SS), where CDW and off-diagonal long-range 
order coexist [11, 12, 14]. 

In the present paper we study the ID BFHM in the limit 
of large fermion hopping which allows for a rather compre- 
hensive understanding of the existing phases and their origin 
in particular in the case of half filling of the spin-polarized 
fermions. In the large hopping limit the fermions can be for- 
mally integrated out [11, 13, 15] resulting in an oscillating 
mean-field potential as well as a long-range density-density 
interaction between the bosons. This interaction has alternat- 
ing sign if the fermion filling is commensurate with the lattice 
which is the origin of the 4k p CDW. In the thermodynamic 
limit it is however formally divergent and needs to be renor- 
malized which is done here by taking into account the back- 
action of the bosons to the fast fermions. The resulting effec- 
tive boson model allows an analytic and quantitative predic- 
tion of the (hb — Jb) phase diagram, where /is is the bosonic 
chemical potential and Jb the corresponding hopping ampli- 
tude. At double-half filling, i.e. Qf = Qb = \, we identify 
an incompressible CDW phase and study its transition to a SF 
with increasing Jb both using analytic results from the effec- 



tive model and numerical simulations based on DMRG [16]. 
DMRG simulations also show the presence of a novel phase 
with coexistence between spatially separated Mott-insulator 
and CDW regions for non-commensurate boson filling. This 
phase which is absent for an extended Bose-Hubbard model 
[17] can be well explained within the effective model and is 
shown to exist for all values of the boson-fermion interaction. 
As the effective theory describes the appearance of a den- 
sity wave on a quantitative level it is expected to explain also 
the conditions for the existence of a SS found in [14] using 
quantum Monte Carlo simulations. Numerical evidence was 
given [14] that a SS exists only if the filling of spin-polarized 
fermions is exactly one half but that of bosons is not commen- 
surate and if the boson-fermion repulsion exceeds that of the 
bosons. The emergence of a finite superfluid fraction from the 
effective model will be discussed at a different place [18]. 

Mixtures of ultracold bosons and spin-polarized fermions 
in optical lattices are well described by the Bose-Fermi- 
Hubbard Hamiltonian [19, 20] 

h = -JbY, + + f £fy (»j - 1) 
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(1) 



where ,b (&, c) are bosonic (fermionic) creation and anni- 
hilation operators and h (m) the corresponding number opera- 
tors. Here, the bosonic (fermionic) hopping amplitude is given 
by Jb (Jf), and U (V) accounts for the intra- (inter-) species 
interaction energy. In the following we consider the limit of 
large fermionic hopping, i.e. we assume Jp S> U, \V\,Jb 
and the energy scale is set by U = 1. 

In this limit of large fermionic hopping the physics of 
the BFHM is well captured by the bosonic phase diagram 
alone. Considering the most interesting case of half filling 
of the spin-polarized fermions, we have plotted in figure 
1 the phase diagram for the bosons obtained numerically 
by DMRG simulations and exact diagonalization (ED) for 
Jp = 10, V = 1.25. Besides the MI and SF phases expected 
from the pure bosonic model, the phase diagram displays a 
CDW phase at double-half filling (gp = qb = 1/2). Exact 
diagonalization is used for very small Jb to avoid boundary 
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FIG. 1: (Color online) Boundaries of the incompressible MI phases 
and the CDW phase for half fermion filling qf = 1/2, Jp — 10, 
and V = 1.25 obtained by DMRG and for small values of Jb by 
ED. One recognizes partial overlap between MI and CDW phase for 
small values of Jb indicating regions of spatial phase separation (PS) 
between MI and CDW. The dashed lines are results from the 2nd 
order perturbation theory based on the effective bosonic model 



effects. Furthermore a novel phase is visible for the case of 
non-commensurate boson filling in which spatially separated 
regions of bosonic Mott insulators and density waves coexist 
(phase separation, PS). Figure 2 shows numerical results 
for the density cuts from within the corresponding phases. 
We note that while the pinning of the CDW to the bound- 
aries is a result of the open boundary conditions required for 
DMRG, the phase separation persits for large systems and was 
verified for small systems using periodic boundary conditions. 
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where g\ = -i^ dr ((Tm J -(r)m J -+;(0)}) F . T denotes 
time-ordering. One recognizes a fermion-induced mean-field 
potential proportional to (rhj)F- The couplings gi describes 
a long-range density-density interaction between the bosons 
separated by I lattice sites. The mean-field potential and 
the density-density interaction are the only interaction terms 
emerging in the effective theory since higher order moments 
in the cumulant expansion vanish. 
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In the following we derive an effective bosonic model 
which provides an understanding of the phase diagram in fig- 
ure 1 on a quantitative level. In the limit of fast fermions one 
could expect that their main influence is through a mean field 
contribution, which according to eq. (1) amounts to a simple 
shift of the bosonic chemical potential \xb — ► Hb — Qf V. 
Indeed the two Mott lobes in fig. 1 are symmetrically located 
around /j, = i V. To explain the CDW and PS phases one 
needs however an effective description beyond the mean-field 
level. In order to adiabatically eliminate the fast fermions we 
first separate the fermion Hamiltonian. As will be seen later 
on it is essential to take into account the back action of bosons 
to the fermions. To do this we incorporate in the fermion 
Hamiltonian the interaction with a mean field potential given 
by a yet undetermined density nj of bosons. Thus = 

-Jf J2j (cjCj+i + + is the fermionic 

Hamiltonian and H\ = V^2j — tne interaction. 

We start from the S-matrix S — T cxp{ — i <1tH[(t)} 
of the full system with Hi[f) = e-^~Hi)r fj^iH-H^r 
and trace out the fermions exactly using a cumulant expansion 
[21]. Since the characteristic time scale of fermionic correla- 
tions is of order 1/Jp and thus much shorter than any other 
time scale in the system, a Markov approximation can be used, 
replacing two-time bosonic operators by equal time operators, 



FIG. 2: (Color online) Densities of bosons and fermions correspond- 
ing to the different regions in the boson phase diagram fig. 1 for a 
lattice of length L — 128 and open boundary conditions. The boson 
number N B is (a) 19, (b) 64, (c) 96, and (d) 101. The plots (a,b,c) 
are taken at J B = 0.01 and (d) at J B = 0.07. While (b) displays 
the gaped CDW, (a) and (c) shows the PS phase, (d) is outside of the 
parameter regime with a CDW. 

In the case of free fermions, i.e. ignoring the back-action of 
bosons (rij = 0), fermionic correlations and couplings gi can 
easily be calculated. For gp = 1/2, gi scales asymptotically 
as 9i ~ (~ 1) /i> i- e - nas a long-range character and alter- 
nating sign. For a general fermion density gp they oscillate 
with period 1/qf which is typical for induced interactions of 
the RKKY type [22-24]. The oscillation of the interaction 
is the origin of the formation of charge density waves. For 
qf = 1/2 the interaction energy is minimized if the bosons 
occupy sites with distance 2. An effective theory with 
coupling constants resulting from free fermions has however 
a fundamental problem: As gi ~ \jl the boson-boson 
interaction energy diverges logarithmically with the total 
length of the lattice. This would result in an incompressible 
CDW for any value of the bosonic hopping Jb- Thus such 
a theory completely fails to describe the transition from a 
CDW phase to a bosonic superfluid. An accurate description 
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of this transition is however important e.g. to explain the 
PS phase and to address existence conditions of a supersolid 
phase. Therefore it is necessary to renormalize the effective 
interaction. This is done here by taking into account the back 
action of the bosons through the mean-field potential in Hp. 
Since for qf = 1/2 the bosonic system is driven into a CDW 
with period 2, a good ansatz is nj = | [l + ( — l) J 7y] , where 
V = v{Jb) is the amplitude of the density oscillation with 
r]\j _ = 1. In general rj is treated as a free parameter and 
can be determined self-consistently by a minimization of the 
energy. For the following calculations it turns out to be more 
convenient to introduce a parameter a proportional to rj as 
a = riV/(4V2nJ F )- 

In order to calculate the fermionic cumulants and therefore 
gi analytically from Hf in the presence of the boson-induced 
mean field potential we make use of Green's function tech- 
niques and perturbation theory in rj V/Jf (for details see [18]). 
This gives gi as a function of a as 
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Similarly one finds for the fermionic density 
1 
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with r]F = 



K 



and K[x] being the complete 



elliptic integral of the first kind. This equation along with 



l + (-l) J 7? 



(5) 



gives the density distributions of fermions and bosons for 
double-half filling, i.e. in the CDW phase, as function of 
the variational parameter a. In the limit a — > the above 
expressions reduces to the free fermion case ((rhj) = h). 
To provide a first test of the validity of the effective theory 
we have plotted in Fig. 3 the ratio of the amplitudes of 
the bosonic and fermionic density waves obtained from the 
data in fig. 1 along with the prediction from eqs. (4) and 
(5). Note that this ratio is exactly fixed by the effective 
theory and independent on the variational parameter a. One 
recognizes a very good agreement. Also shown in Fig. 3 
are the amplitudes of bosonic and fermionic density waves 
respectively as function of bosonic hopping Jb obtained 
numerically as well as from the effective theory using the 
numerical data of the other species as input. It can be seen, 
that the renormalized effective theory fits quite well with the 
numerical results. 

Using the effective Hamiltonian, (2) we will now derive an 
analytic approximation to the phase diagram of the full BFHM 
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FIG. 3: (Color online) CDW amplitudes (points) of bosons and 
fermions for double half filling as function of normalized boson hop- 
ping. Dash-dotted line: calculated factor a from the fermionic data. 
Dashed line: analytic estimate of the bosonic amplitude from the 
fermionic input a. One recognizes a very good agreement between 
the numerical data and the analytic estimate, where it should be kept 
in mind that the underlying perturbation theory gets better for rj — > 0, 
i.e. for increasing Jb ■ 



using a strong-coupling expansion valid for small values of 
Jb [26]. Since we are mainly interested in the boundaries 
of the incompressible lobes, we will calculate them in the 
canonical ensemble from the energies of the relevant states as 
a function of the bosonic hopping amplitude Jb- The upper 
(lower) boundary is given by the bosonic particle-hole gap. 
With this, the chemical potentials for bosonic filling qb are 
given by M ± B = ± (E(g B L ± 1) - E{g B L)), where E(N) 
is the ground state energy for a given number of bosons N, At 
Jb = 0, this is straight forward to calculate since the ground 
state distribution of the bosons is trivial and the variation pa- 
rameter 77 is fixed to unity. 

For Jb > we apply degenerate perturbation theory in 
2nd order using Kato's expansion [25]. In second order, there 
is a local correction of the ground state energy for all numbers 
of particles, as well as 2nd order two-site hopping processes 
connecting the states within the ground state manifold in the 
case of an additional (absent) boson. Incorporating this, the 
upper and lower chemical potentials for the CDW phase can 
be expressed in a simple analytic form as 



, V 

Mi = tt ±Vr)F±go(a) - (3±J B . 
2 2 



(6) 



Similarly one finds for the chemical potential corresponding 
to unity and zero filling 



V 

Mi = ~2 - -9o(0) + 2J B - aJ B , 
(4 = ^+.9o(0)-2J B . 



(7) 
(8) 



Here rjp and go{a) are taken for rj = 1 (Jb = 0) because 
of the perturbation expansion. The derivation of a and j3± is 
lengthy but straight forward. Their explicit form will be given 
along with other details elsewhere [18]. One finds in partic- 
ular (3 + > > Note that Vr\p is positive irrespective 



4 




interaction V 

FIG. 4: (Color online) Boundaries (shifted by V/2) of incompress- 
ible MI and CDW phases for half f ermion filling qf = 1 /2 at vanish- 
ing bosonic hopping Jb = with varying interaction V. Curves are 
the analytic results and the data points are obtained by DMRG (CDW 
phase, L = 128) and ED (MI lobes, periodic boundaries with infi- 
nite size scaling). The shaded region depicts the coexistence phase 
between MI and CDW with fig > /i7 and /xj > fa . 

2 2 

of the sign of V and is larger in magnitude than both go(0) 
and go{a). With this one can see that /it I _ n > Hi I _ n > 

u,7\ , n < Un\ , „• Thus there exists a region where 

the chemical potential is not monotonous in the boson num- 
ber. This explains the coexistence of MI and CDW in the PS 
phases found numerically in fig. 1. The long-range character 
of the fermion mediated interaction together with the fermion- 
induced mean-field potential prefers extended, spatially ho- 
mogeneous regions of a commensurate CDW. Extra bosons 
will be pushed out and form an incompressible Mott insulator 
region. Such a phase does not exist in purely bosonic systems 
with extended interactions due to the absence of the fermion 
induced mean-field potential V t\f- If the bosonic hopping 
exceeds a certain critical value, given by the crossing of the 
curves /it with /t^ or respectively /i^ the minimization of ki- 

2 

netic energy by equally distributing the particle is larger than 
the loss in interaction energy due to the fermion mediated in- 
teraction. 

In fig. 4 we have plotted the chemical potentials for zero 
bosonic hopping as function of the interaction strength V 
defining the boundaries of the Mott insulating phases with 
zero and unity filling as well as the lower and upper bound- 
aries of the CDW phase with half filling of bosons. One rec- 



ognizes that phase separation between MI and CDW exists 
for all values of the boson-fermion interaction V. Once again 
there is a rather good agreement between full numerics and 
effective theory, which provides another test for its validity. 

The phase boundaries for Jb > obtained from the 
analytic results for the chemical potentials in eqns. (6) to (8) 
are shown in fig. 1 as dashed lines. Although the precise form 
of the CDW lobe is not correctly reproduced (as expected for 
the strong-coupling perturbation approach), the qualitative 
agreement is remarkable. The strong coupling approximation 
yields a critical value of J^ DW = 0.025 beyond which the 
CDW ceases to exist for J F = 10 and V = 1.25. Whether or 
not the CDW gap vanishes at a finite value of Jb is however 
unclear. Our numerics indicates that a very small gap may 
persist even for values of Jb = 1 and may not close at for 
Jb < Jf- The critical values Jg S w 0.01 for the PS region 
obtained from the effective model agrees however rather well 
with the numerical data. 

In summary we developed an effective model for a mixture 
of bosons and spin polarized fermions in a periodic lattice in 
the limit of large fermion hopping. This model reveals the 
physical origin of the incompressible CDW phase and pro- 
vides a simple quantitative description. The fast fermions 
cause a mean-field potential and mediate a long-range density- 
density interaction which is of alternating sign for gp = 1/2. 
In order to accurately describe the conditions for the existence 
of a CDW renormalization effects due to the back-action of 
the bosons need to be taken into account. The density wave 
amplitudes where obtained from an analytic model and ver- 
ified by numerical DMRG simulations. The effective model 
also gives a simple understanding and quantitative description 
of a novel phase where spatially separated regions of a maxi- 
mum amplitude CDW and a MI coexist. The effective model 
is expected to provide a means for predicting and understand- 
ing conditions for the existence of a SS phase in Bose-Fermi 
mixtures and other mass imbalanced two-species models. 
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